arXiv:1505.03004v3 [cond-mat.str-el] 11 Apr 2016 


Light induced magnetization in a spin S = 1 easy-plane antiferromagnetic chain 
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The time evolution of magnetization induced by circularly polarized light in a S' = 1 Heisenberg chain with 
large easy-plane anisotropy is studied numerically and analytically. Results at constant light frequency D = Do 
are interpreted in terms of absorption lines of the electronic spin resonance spectrum. The application of time 
dependent frequency D = D(f) light, so called chirping, is shown to be an efficient procedure in order to obtain 
within a short time a large, controlled value of the magnetization . Furthermore, comparison with a 2-level 
model provides a qualitative understanding of the induced magnetization process. 
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Far-from-equilibrium condensed matter physics is a chal¬ 
lenging, still largely uncharted territory. Concerning the out 
of equilibrium dynamics of quantum magnets, the control 
of magnetic properties by means other than a conventional 
magnetic field is of strong current interestFor instance, 
engineering the quantum state, i.e. wavefunction, is essen¬ 
tial for quantum simulators, precision sensors or spintronic 
devices^*®. Recent experimental advances allow to manipu¬ 
late the elementary low-energy excitations with terahertz laser 
pulses" a prominent example being the ultrafast coherent 
control of antiferromagnetic magnons. A time-dependent ( 
rotating) magnetic field of highly intense terahertz laser pules, 
with photon energy below the electron energy scale, con¬ 
trolled the coherent spin waves without interfering with the 
motion of charge carriers. 

In quantum magnets with reduced dimensionality the ther¬ 
modynamic and transport properties exhibit a rich magnetic 
field dependence related to the total magnetization of the 
system. Prominent examples of such a behaviour are the field- 
induced quantum phase transitions of the organic compound 
NiCl2-4SC(NH2)2 (dichlorotetrakisthiourea-nickel abbrevi¬ 
ated as DTN). At zero temperature, the first transition oc¬ 
curs at a critical field hi where the energy gap closes and 
a finite magnetization develops in the ground state (GS); 
the second one occurs at /12 where the magnetization fully 
saturates leading to a ferromagnetic GS. By now, the low- 
energy physics of the DTN compound has been well studied 

experimentally '^ 42-26 understood theoretically. The ba¬ 
sic model that describes the magnetic excitation spectrum of 
DTN was found to be the one dimensional S' = 1 antiferro¬ 
magnetic Heisenberg model (AHM) with exchange coupling 
constant J and large easy-plane anisotropy D. As shown 
in Refs. 27-30 such a Hamiltonian reproduces in great detail 
the low lying electronic spin resonance (ESR) spectrum. The 
anisotropy D/ J ^ 4 of DTN, being the largest energy scale 
in the system, is responsible for a large energy gap 0{D) that 
can be closed by a magnetic field h. 

In this work we study the rotating magnetic field induced 
nonequilibrium magnetization in large, easy - plane 
anisotropy AHM. For a field rotating at constant frequency 
(circularly polarized light) we are connecting the numerical 


results with the linear response (LR) theory predictions for the 
transition frequency of the corresponding ESR experiment. In 
the case of a chirped (time dependent) frequency of the light^', 
our results indicate that the short-time behaviour of the mag¬ 
netization is mainly driven by the anisotropy part of the sys¬ 
tem. This time scale, together with the dependence of the 
magnetization on the chirp parameters, can be accurately de¬ 
scribed by a 2- level model. The dynamics beyond the charac¬ 
teristic time of the latter is dominated by the Heisenberg part 
of the model. Although we focus on DTN as a typical one di¬ 
mensional S' = 1 easy- plane AHM, our analysis is also valid 
for other Hamiltonians, e.g. a 2-level model will yield the 
correct physics for S = 1 models in all dimensions provided 
that D ^ J. 

As prototype model we choose the S = 1 AHM with 
single-site, easy-plane anisotropy Z? on a chain with L sites 

L 

Ho = Y, • S*+i + D{S!)^ + hS!] , (1) 

i=l 

where = (Sf, Sf, Sf) are spin S = 1 operators at site i, 
Sl +1 = Si (periodic boundary conditions), h is a magnetic 
field, and J(^ 2K) the antiferromagnetic exchange constant 
(we will further on use h = ks = Pb = 1 and set J = 1 as 
the unit of energy). Hereafter, we will use £> = 4 8K) and 

for such an anisotropy the critical fields are: hi ~ 2.28 and 
/i 2 = 8^^. We will assume that only the magnetic component 
of light, propagating in the z-direction, couples to the system. 
The time- dependent Hamiltonian of the corresponding setup 
can be written as 

L 

Hit) =Ho-AJ2 , (2) 

i=l 

where A > 0 and H > 0 are the amplitude and frequency 
of light respectively and Sf are spin raising and lowering op¬ 
erators. Thus, each spin “feels” a magnetic field rotating in 
the xy-plane, 2A ^JS'f cos(Hf) + S'f sin(Hf)]. The mag¬ 
netization induced is positive, in order to obtain a negative 
magnetization one should substitute H —Cl in (2). Note 
that in a real experiment a propagating light pulse has some 
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time and frequency dependence, an issue that we will discuss 
later on. In order to probe the sample magnetization perpen¬ 
dicular to the polarization plane one can use a second optical 
pulse and measure the change in its polarization state induced 
by the magnetization either in transmission (Faraday effect) 
or reflection (Kerr effect) geometry. 

The time evolution of the magnetization is given by 


M\t) 


(vl/(f)|54|vh(f)) 


(3) 


where = {^/L)Y^-S^ and is a solution 

of the time-dependent Schrodinger equation 6c)t|T'(f)) = 
H(t)). In our calculation we choose 6t in such a way 
that typically (T'(f)l'I'(f)) ~ 1 at any time t {6t ~ 10“^). A 
general procedure goes as follows: (i) first, with help of exact 
diagonalization we calculate the GS of(l), |'I'(—5f)) = |GS). 
(ii) Next, at time f = 0 we instantaneously turn on the light 
and (iii) finally, we perform the time evolution of it on the ba¬ 
sis of the time-discretized version of the Schrodinger equation 
with (2) (using a fourth-order Runge-Kutta routine). 

Let us first focus on the system (2) at constant frequency 
n = rig. It is clear that the maximum value of the magne¬ 
tization is induced by light at the resonance frequency of the 
system, fig = ^r, that can be interpreted in the spirit of an 
ESR spectrum. For small enough A (<C J) the system is in 
the linear response regime and the low absorption lines of the 
ESR spectrum of (l)^®’^^ correspond to the resonance frequen¬ 
cies riR of (2) at given h. Eurthermore, (2) at fl = fig can be 
mapped by a unitary transformation (or Eloquet theory) to an 
effective static model^'*^^®, where the latter has a form similar 
to the one when dealing with an ESR experiment. Note that 
the same procedure was used in Ref. 34 in order to study a sys¬ 
tem with small magnetic anisotropy D = 0.25 (Haldane-like 
limit). 

Eig. 1(a) depicts a typical example of the time dependence 
of as a function of time for a system with h = 0 and 
constant H. Several conclusions can be drawn directly from 
the obtained results: (i) It is evident that the induced by 
Hg = is dominating above other frequencies, (ii) The 
beating frequency presented in the inset of Eig. 1(a) is at¬ 
tributed to finite size effects, (iii) The value of D,r = 6 for 

= 0 is consistent with the lowest transition lines of ESR 
spectrum. In fact, in the gapless regime h < hi, the ESR 
lines can be calculated by a 1 /D expansion^^’^^, i.e.. 


ujji — D h, ujR — D -\- ‘2J — h . 


Such lines correspond to transitions from the GS to states with 

= ± 1 . 

In Eig. 1(b) we present a heat map of the average (over time 
span dt < t < 100) net magnetization, — M^{t = 0), as 
a function of magnetic field h and frequency fig. Our results 
perfectly reproduce both ESR predictions, e.g. see Eig. 6 of 
Ref.^^. In the considered held h region we also see continu¬ 
ation of the OJG = D + h line - transitions from a magnon 
to a single-ion bound state. Other resonance lines can also be 
captured, e.g. transitions from the fully ordered ferromagnetic 
state in the h > h 2 region, can be resolved by looking for Q,r 



Figure 1. (Color online) (a) Magnetization as a function of time 
M^{t) calculated for L = 11, A — 0.1, h = 0, and fig = 4, 6, 8. 
Dashed horizontal line represents average value for fig = 6. Note 
that the results for fig = 4 and 8 are multiplied by factor of 5 for 
clarity. Inset: M^{t) induced by fig = Q,r = 6 (as in the main 
panel) for f up to f = 250. (b) Heat map of average net magne¬ 
tization, = 0), as a function of magnetic field h and 

frequency fig, calculated for L = 10 , A = 0.1. ujb line (red color 
in the heat map) is obtained with fl > 0 in (2), uja,g (blue color) 
with fl < 0. Solid and dashed lines represent the liJa,b,g ESR res¬ 
onance lines and their continuation into the gapless regime. Vertical 
solid line represents the critical field hi. 


of negative magnetization. In Eig. 2(a) we present = 5) 
as a function of frequency Hg. The maximum value of mag¬ 
netization for given h and H > 0, or H < 0, is consistent with 
the ESR predictions. 

Although we chose the GS as the starting point of the time 
evolution this is not a zero temperature (T = 0) result. Within 
LR theory one would expect for T = 0 rather sharp transition 
lines^®. It is clear from Eig. 1(b) that our resonance lines are 
not (5-peaks, with nonzero intensity for all considered transi¬ 
tions oja,b,g- Also, in Eig. 2(b) we present the dependence of 
h = 0 average magnetization on frequency fig for various 
amplitudes A = 0.01,0.05,0.1,0.5 in (2). Within LR such a 
broadening of the line could be interpreted as the increase of 
an effective temperature. 

Next, in order to induce a macroscopic magnetization in a 
controlled way we study the application of a chirped pulse, 
fl = n{t). Although the time dependence of fl can be com¬ 
plicated and its functional form dependent on the experimen¬ 
tal setup, the main features should be captured by the simple 
form 


fi(<) = rii — vt, 


(4) 
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Figure 2. (Color online) (a) Frequency flo dependence of the mag¬ 
netization at time t = 5 and L = 11. Results for h — 0, 2,4 are 
calculated with flo > 0 (positive magnetization) and for h — 2,8 
with flo < 0 (negative magnetization). Note that for h > hi (pre¬ 
sented for h = 4,8) the ground state has net magnetization already 
at f = 0. (h) Frequency dependence of average magnetization 
for h = 0, L = 11 and various amplitudes A = 0.01, 0.05, 0.1, 0.5. 
Results for A = 0.01 are multiplied hy factor of 5 for clarity. 
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where £7/ is the initial (t = 0) frequency and v is the chirp, 
i.e., the “speed” of frequency change. Within such a notation 
ilj = ri/j and 1 / = 0 corresponds to a time independent 17 at 
the resonance frequency. In the following, we will consider 
only the h = 0 case, i.e, iln = 6, as we would like to study 
the magnetization induced only by light. 

The qualitative dependence of the magnetization on the am¬ 
plitude A and chirping iz can be understood within a 2-level 
model. 


i ?2 =0|0)(0|-F77|l)(l| + y2A(e-^^*|l)(0|+h.c.) , (5) 


where |0) (|1)) correspond to the = 0(1) states of the term 
D{Sf)'^, relevant to the J / D —?> 0 limit of (1). Note that the 
resonance frequency of (5) is simply VIr = D. Within this 
model, a perturbative expression, a = Aj^ 0, of the 

time dependence of the magnetization can be given as. 


l + |FF(t)|2’ 


I 

Wit) = V2A J 
0 

t^/y 


( 6 ) 


Figure 3. (Color online) Time dependence of the magnetization cal- 
culatedforL = 11, h = 0, A = 0.1. (a) Magnetization as a function 
of time calculated for various initial frequencies 17/ = 6,8,10 and 
v — 0.01. The horizontal line represents the average value of mag¬ 
netization for v — 0 and 17/ = 17// = 6 (resonance frequency), (b) 
Comparison of the magnetization as calculated with Eq. (3) for, (i) 
the full Hamiltonian (2) with 17/ = 8, iz = 0.01, (ii) the 2-level 
model (5) and the perturbative solution - Eq. (6) with A = 2, 
iz = 0.01. (c) Magnetization as a function of normalized time vt for 
iz = 0.005, 0.01, 0.05, initial frequency 17/ = 8 and A = 0.1. Inset: 
the same results as a function of time t. 


where A = 17/ — Hr. It is obvious from the above equation 
that M^{t) depends only on a and the detuning A. 

In Fig. 3(a) we present the magnetization dependence on 
the initial frequency 17/ at fixed iz. In panel (b) we present 
numerical results obtained from the full Hamiltonain (2), the 
2-level model (5) with corresponding detuning, together with 
the perturbative solution Eq. (6), that captures the main fea¬ 
tures of the magnetization profile. Note that the results are 
indistinguishable till the saddle point of Eq. (6), i.e., at tg = 
A/2zz. Erom the results presented in Eig. 3(b) it is obvious 
that the main effect of the exchange coupling J is in the dy¬ 
namics of the magnetization at times t beyond tg. It is also 
interesting to note that the magnetization induced by a con¬ 
stant frequency light 17 = Hr, as indicated by a dashed line 
in Eig. 3, is reached at the saddle point time tg. We observe 
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such a behaviour for all flj > flu. 

In Fig. 3(c) we show the time dependence of the induced 
magnetization for different chirping speeds ly. We observe 
that, (i) in the scaled time vt the curves are practically identi¬ 
cal with the crossing of the mean value at = 0 and flj = flu 
at lyt ^ vts = 1, reaching maximum at vt ~ 2; (ii) the 
magnetization at long times is weakly dependent on time, a 
remarkable result considering that we are dealing with a full 
many-body problem, where a decay could be expected; (iii) 
it is clear that, as the total magnetization com¬ 

mutes with the Hamiltonian, after switching off the light at a 
certain time, remains constant at its instantaneous value. 
This allows for a tight control of the value of the induced 
magnetization in the system. Further simulations for differ¬ 
ent fli > flu confirm this picture; crossing the resonance 
frequency by chirping the light frequency induces a stable 
macroscopic magnetization in the system. Additionally, it 
is clear from the solution of the 2-level model that inverting 
A' = —A and ly' = —v produces identical evolution of the 
magnetization. 

Considering the amplitude and chirping speed dependence 
of the long time asymptotic magnetization achieved, first of all 
we observe that the 2-level model can be mapped in a rotating 
frame to a Landau-Zener type tunneling problem, 

H 2 = A|T)(T| - A|0)(0| + (|T)(0| + h.c.) , (7) 

where A = A/2 — In the Landau-Zener problem A = 0 
and the time evolution is from t = —oo to t = -foo, while 
in the situation we are considering the time evolution starts 
at f = 0 and from a finite frequency shift A. For A/iy ^ 
1 the asymptotic M^{oo) coincides with the probability of 
occupation of level |1) given by the Landau-Zener expression 
1 — exp(—7ra^). 

In Fig. 4 we present a comparison of the long time mag¬ 
netization {i/t = 2) in the full model (2), the perturbative pre¬ 
diction Eq. (6) and the Landau-Zener expression. Note that al¬ 
though is a perturbative solution {a 0) and the detailed 
dynamics beyond tg is not captured correctly (see Fig. 3(b)), 
the overall agreement of the asymptotic magnetization is qual¬ 
itatively described till a ~ 136. 

Finally, turning to the experimental realization, for the DTN 
compound (J = 2.2 K, D = 8.9 K) the resonance frequency 
is flfi Ri 300 GHz. Light of magnetic field intensity Ri 
0.3 Tesla corresponding to an electric field of Ri 1 MV/cm and 
a chirping speed i/ Ki 0.1 will induce a controlled macroscopic 
magnetization within Ri Ipsec^^. In a realistic experimen¬ 
tal situation, several issues arise, (i) in terahertz spectroscopy 
the light is in the form of a pulse of duration Ri 1 psec, (ii) 
the effect of electric field should be estimated, (iii) experi¬ 
ments are at a finite temperature, (iv) there is spin-lattice re¬ 
laxation which could be detrimental to the process of inducing 
a macroscopic magnetization. However, it is known that in 
several quantum magnets^® the relaxation time is surprisingly 


long. Preliminary finite temperature simulations and consider¬ 
ations are encouraging in rendering the proposed experiment 
feasible. We should also note that the large variety of quantum 



Figure 4. (Color online) Scaling parameter a = Aj^/v dependence 
of the magnetization as calculated for L = 11, i/ = 0.1 (full squares) 
and 1 / = 0.01 (open squares). The snapshot of magnetization of the 
full model (squares) is taken at i/t = 2, i.e, f = 20 for i/ = 0.1 and 
t = 200 for v = 0.01. Circles represent the magnetization M^(oo) 
dependence on a in the 2-level model (5) and the black dashed line 
depicts the Landau-Zener expression. 


magnets, allow for a tailoring of the experiments in terms of 
light frequency, relaxation time etc. 

In summary, we have studied an efficient protocol which in¬ 
duces magnetization without external magnetic field applied 
to the system. Results for circularly polarized light pulse 
at constant frequency are explained with the help of reso¬ 
nance lines of ESR transitions at finite temperature. We have 
also presented comprehensive results on the dependence of 
the magnetization on a chirped pulse. The latter, experimen¬ 
tally relevant, protocol can be qualitatively and even for some 
time scales quantitatively described with the help of a 2-level 
model. Also, it was shown^® that (1) can be mapped to an ef¬ 
fective S' = 1/2 AHM with exchange anisotropy < 1. Our 
2-level predictions for this model will be even more accurate 
since the mapping favors the large-!? limit. 
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